A flexible method to calculate the distributions of discretization errors in operator-split codes with stochastic noise in problem data

This report investigates a technique to calculate the distributions of discretization errors for a model of advection-diffusion-reaction with stochastic noise in problem data. The focus is on operator-split discretization methods. The error is decomposed into components due to the splitting and due to the discretization within each component. We present a method to estimate the distributions of the total discretization error for a quantity of interest as well as the aforementioned error components. These computations require information about the propagation of uncertainty within and between code modules, for which purpose we use a flexible code framework for uncertainty quantification. Computational examples illustrate that in certain cases the distribution of numerical error can have a significant impact on the code output. Using information about the components of the numerical error, decisions can be made regarding discretization parameters to efficiently reduce the impact of numerical error on the distribution of model outputs.